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Quantum tunneling through a two-dimensional static barrier becomes unusual when a momentum 
of an electron has a tangent component with respect to a border of the prebarrier region. If the 
barrier is not homogeneous in the direction perpendicular to tunneling a fraction of the electron 
state is waves propagating away from the barrier. When the tangent momentum is zero a mutual 
interference of the waves results in an exponentially small outgoing flux. The finite tangent momen- 
tum destroys the interference due to formation of caustics by the waves. As a result, a significant 
fraction of the prebarrier density is carried away from the barrier providing a not exponentially 
small penetration even through an almost classical barrier. The total electron energy is well below 
the barrier. 



PACS numbers: 03.65.Sq, 03.65.Xp 



I. INTRODUCTION 

The problem of tunneling through a potential barrier 
was a subject of attention for many years in nuclear 
physics, particle physics, atomic physics, and chemical 
physics. Quantum tunneling across a one-dimensional 
static potential barrier is described by the theory of 
Wentzel, Kramers, and Brillouin (WKB) Q if the barrier 
is not very transparent. Early studies of quantum tun- 
neling were performed in Refs. 0, S Some modern 
aspects of thephenomcnon can be found, for example, in 
Refs. |,i,0,i. 

A scenario of tunneling through a multi-dimensional 
barrier is well studied @, [lol, lill, [H, [S El E d, 
[itI . [TsI . [19I . |20| . The main contribution to a tunneling 
probability comes from the extreme path in the {x, y} 
plane (we consider two dimensions) linking two classi- 
cally allowed regions as drawn in Fig. [Tfa). The path is 
a classical trajectory with real coordinates which can be 
parameterized by imaginary time. The underbarrier tra- 
jectory is a solution of Newton's equation in imaginary 
time. The trajectory is given rise by a particle hitting 
with a zero tangent momentum a border of the classi- 
cally allowed region from the the classical side. This is 
shown in Fig. [ija). Under the barrier the probability 
density reaches a maximal value at each point of the tra- 
jectory along the orthogonal direction with respect to it. 
Therefore around the trajectory, which plays a role of a 
saddle point, quantum fluctuations are weak. The wave 
function, tracked along that trajectory under the barrier, 
exhibits an exponential decay generic with WKB behav- 
ior. 

However, in some cases tunneling through multi- 
dimensional barriers occurs according to a different sce- 
nario which is far from being similar to WKB. In this 
paper unusual tunneling through two-dimensional barri- 
ers is studied. 

When the prebarrier state has a tangent component 
of a momentum, as in Fig. Hlb), there are no extreme 
points at the border of the prebarrier region since the 



derivative of a wave function along the border is finite. 
This means that a tunneling probability is no more de- 
termined by the main underbarrier path but comes from 
a wide set of paths indicated in Fig. [T]^b) by the dashed 
arrows. Traditionally, a decay of a state with a tangent 
momentum is not considered since it does not correspond 
to a saddle point and, hence, the net contribution is sup- 
posed to be averaged down to a small value. As shown in 
the paper, that conclusion is not correct and states with 
tangent momentum can play a crucial role in tunneling 
processes. 

The phenomenon can be explained in terms of the fol- 
lowing arguments. 

When a two-dimensional barrier is homogeneous per- 
pendicular to the tunneling direction x (Vo(a:)) tunneling 
is generic with WKB mechanism. Suppose that there is 
an impurity localized at the barrier region and described 
by the potential u(x, y). The underbarrier wave function 
in the total potential Vo(a;) + u{x,y) contains a set of 
overbarrier propagating waves if to treat them as eigen- 
functions of Vq{x). When a tangent momentum is zero, 
as in Fig. (Ha), the energy distribution of the propagating 
waves is smooth leading to a strong mutual cancellation 
of their contributions to a total wave function due to in- 
terference. As a result, a contribution of the propagating 
waves to an outgoing flux is exponentially small. 

In contrast, when a tangent momentum is finite, as in 
Fig. [H^b) , the mutual cancellation of those waves due to 
interference can be not complete. Analogous to optics, de 
Broglie waves may form caustics where a distribution of 
waves becomes not smooth [iM [3 IHI ■ This violates 
the strong mutual cancellation of the propagating waves 
and the surviving fraction of them goes away from the 
barrier providing a not small output. More details are 
given in Sec. IXIII 

The underbarrier interference opens a possibility of 
penetration through almost classical barriers. The en- 
hanced tunneling through a static barrier was initially 
proposed in Ref. [2^ . This phenomenon was also stud- 
ied in tunneling through nonstationary barriers [2^ . [zj . 
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|25| . l26l . |27| where an underbarrier phase was created by 
quanta emission. 

An exact wave function in tunneling problem can be 
written in the form C{x,y)exp[iS{x,y)/h] where 5 is a 
classical action and C is a prefactor. An exact analytical 
solution of the Schrodinger equation for a complicated 
barrier is impossible. We use a semiclassical approach 
when only the exponential part of the wave function is 
calculated. This is called the main exponential approx- 
imation when the action S is large and the prefactor C 
is less important. The classical action is calculated in 
Sec. I VII as a solution of the Hamilton-Jacobi equation. 
The action was also tracked along classical trajectories 
in imaginary time. The both methods lead to the same 
tunneling probability. 

At the first sight, the obtained solution of the 
Hamilton-Jacobi equation provides the main (exponen- 
tial) part of the exact wave function with no problems. 
But there is a delicate phenomenon in quantum mechan- 
ics which may destroy a semiclassical solution. When 
two branches of a wave function exist at the same spa- 
tial domain the smaller one can unexpectedly disappear. 
From the semiclassical standpoint it happens as a jump 
in space of the prefactor C(x. y) d own to zero. This is 
called Stokes phenomenon [281. l29l[30t . In Sec. IVIII it is 
shown that there are no Stokes jumps in the semiclas- 
sical solution obtained. So despite the exact solution is 
unknown, it is proved that exp[z5(a;, is the main 

part of the exact wave function and the prefactor plays 
a secondary role. 

In principle, a numerical solution of the static 
Schrodinger equation in two dimensions is possible. But 
in our case one should keep at some spatial domains two 
branches when one of them is exponentially small. In or- 
der to resolve that branch a precision of calculations has 
to be very high resulting in an extremely large number 
of discrete points in a two-dimensional net using for a 
numerical computation. This essentially exceeds mem- 
ory facilities of computers. So a numerical study of the 
problem seems to be impossible at present. 

In Sec. IXII tunneling from a quantum wire through a 
barrier with an impurity is investigated. An influence 
of impurities inside a potential barrier on tunneling was 
widely studied. See, for example, [U, [H, [H, [s^l • A fa- 
mous mechanism is resonant (Wigner) tunneling when an 
impurity level coincides with a particle energy This 
is not our case since we are away of Wigner resonance. 
Another famous mechanism is the interference in scat- 
tering in a system of many impurities as in localization 
phenomena |35|. This is also not our case since there is 
just one underbarrier scattering center. 

The proposed mechanism of impurity assisted tunnel- 
ing was in shadow in the previous studies. A neces- 
sary condition is a tangent momentum in a prebarrier 
region. Underbarrier waves with the tangent momen- 
tum are scattered by the impurity resulting in overbarrier 
propagating waves which carry away a significant fraction 
of the prebarrier density. 







FIG. 1: Solid curves correspond to a constant potential en- 
ergy. Tunneling occurs between two classically allowed re- 
gions, the prebarrier one is to the left, (a) The conventional 
mechanism. The main trajectory in the classically forbidden 
region is indicated by the dashed path, (b) The electron hits 
the border with a tangent velocity. The set of subsequent 
propagating waves is shown by the dashed arrows. 



II. FORMULATION OF THE PROBLEM 

Below we consider a simpler barrier than one in Fig. [T] 
Tunneling occurs in the x direction from a straight in- 
finitely long quantum wire aligned along the y axis. A 
particle, localized in the {x, y} plane, is described by the 
Schrodinger equation 



2m \ dx"^ 



Vl + /3'(2/)'5(x)V 
~£\x\i! = Ei;. (1) 

The S well in Eq. ([1]) relates to the quantum wire placed 
at the position a; = 0. When the wire is homogeneous 
{/3{y) — 0) the discrete energy level in the S well is 
—uq. The non- homogeneity of the wire is supposed to 
be localized at a finite region of y where /3(y) is not 
zero and therefore /3(±oo) = 0. We also suppose that 
f3{y) = f3{—y). Tunneling occurs through the triangu- 
lar potential barrier created by the static electric field £. 
The energy is below the barrier (E < 0) and in order to 
move from the wire (x — 0) to an infinite x one should 
pass through the potential barrier. For simplicity we use 
the even potential — f to get the problem symmetric 
with respect to x. One can consider positive x only if to 
put the boundary condition 



dip{x,y) 



dx 



^/2rmIo 



^i + PHy)my), (2) 



x=0 " 
which follows from Eq. ([T]). 



III. ZERO ELECTRIC FIELD 

First, let us find a wave function for a wide barrier 
when the electric field £ is zero. In this case the barrier 
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becomes rectangular. By Fourier components the 
function reads 



wave 



dk 

ZTT 



ik{y-yi) - x\ k 



2mE 



(3) 



An equation for '^'{Q, y) can be obtained in the same man- 
ner through Fourier components from Eq. ([2]) 



j dyi^{0,yi) J 



27T ^ 



^ (4) 



= v/i + /3'(y)V'(o,y). 



Eqs. Q and 111 provide an exact description of exponen- 
tial decay under the rectangular barrier. 

In the semiclassical limit the solution of Eq. (j4|) 



'0(0, y) ~ exp 



^^I'dy^dl 



E 

Uo 



(5) 



is obtained with exponential accuracy by a saddle 
method. Validity of the solution ([5]) can be easily checked 
if to insert it into Eq. ^ and to take saddles subsequently 
with respect to yi and k. 

The physical meaning of a state in the non- 
homogeneous wire when an electric field is zero becomes 
clear in the limit when /3(j/) ^ 1 and (E + uq) <C uq- In 
this case the wave vector k is small and Eq. ^ reads 



'4!{x,y) = '0(0, y) exp ("^V^mwo 
Eq. (H)) becomes of a Schrodinger form 

d^^j{Q,y) 



(6) 



2m dy^ 



-«o/3'(y)V(0,y) = (ii; + 7.o)V'(0,z/). (7) 



When E > —uq a, solution of Eq. d?]) at a large negative 
y is a plane wave which reflects from the effective poten- 
tial —uoPHy)- This case is generic with an overbarrier 
reflection At |y| ^ oo 



■0(0,2;) = exp 
+9{-y)y/Rexp 



'-l^2m{E 
n 

iy 



Uo) 



y/2m{E + uo) 



(8) 



When the typical distance y ^ a, where P{y) is not zero, 
is sufficiently large the reflection coefficient \R\ in Eq. ([5]) 
is exponentially small [l| 



\R\ ~ exp 



-c-\/2m{E + uq) 



< c- 1. 



(9) 



IV. FINITE ELECTRIC FIELD 

Below we use the dimensionless form of the 
Schrodinger equation ([1]) measuring x and y in the units 
of uq/£ 



1 fd^ijj d'^i) 



B 



^l + a-'{y)5{x)^ (10) 
-\x\'il) = (7- l)'ip. 



The large semiclassical parameter is 



B = 



uoy/2muQ 



(11) 



a{y) = f3{yuo/£), and the energy is E — (7— 1)1*0. Since 
the energy E is below the barrier the condition 7 < 1 
holds. The equation pH)) can be considered at a; > 
if to impose the boundary condition ^ which in the 
dimensionless units reads 

dtp{x,y) 



-B^l + a^y)^{0,y). (12) 

At positive x the variables in Eq. (jlOp are separated and 
a solution can be written in the general form 

^{x,y)= / dke:^Y>[iBF{k)] (13) 
Jc 

exp (^By^J^ + k'^ - B j dxi\/\ + B ~ x^ . 

The integration contour(s) C lies in the plane of com- 
plex variable k. This relates to Laplace's method for 
differential equations. The function F{k) should be de- 
termined from the condition (I12p . The resulting equation 
for -0(0, y) — j dk exp{iBF+iByy^ 7 + k^) is generic with 
Eq. Q. 

A calculated function F{k) has to be inserted into 
Eq. I|13p. Since the parameter B is large one can use 
a saddle point method to evaluate the k integral in 
Eq. For each x and y this method gives a cer- 

tain saddle value ks{x,y). Quantum fluctuations around 
the saddle determine a preexponential factor. As shown 
below, the saddle formalism corresponds to a general in- 
tegral of Hamilton- Jacobi equation considered in Sec.lVl 



V. HAMILTON-JACOBI FORMALISM 

In this section we develop a semiclassical method of 
solution of Eq. (|10p which is reduced to the Hamilton- 
Jacobi formalism. The wave function is of the form 



ipix,y) = exp [iBa{x,y)] . 



(14) 



If to insert the expression into Eq. PH)) we obtain at 
positive X the equation 



dx 



dy 



B ydx"^ dy"^ 



(15) 



4 



Since i? is a large parameter Eq. is reduced to 



dx 



dy 



7- 



(16) 



Eq. p6)) is called the equation of Hamilton- Jacobi for the 
classical action S = hBa{x,y). The function S{x,y) is 
really classical one since it does not depend on Planck's 
constant as B is inversely proportional to h. Eq. (fTB]) 
holds at positive x. The boundary condition follows from 
Eq. (E 



da{x,y) 



dx 



ly^l + a^y). 



(17) 



x=0 



A general integral of the Hamilton- Jacobi equation can 
be obtained by the method of variation of constants [s^ . 
The general integral, satisfying the condition pT|) . has 
the form 



a{x,y)=i dxi\/a'^[iv{x,y)] + 1 



Xl 



+yV"f + a'^[iv{x,y)] - 



2/1- 



dyi 



where the function v{x,y) obeys the equation 
v{x,y) + iy /"^ dxi 



\/j + a'^[iv{x,y)] Jq y^a'^[iv{x, y)] -|- 1 - xi 
In this formalism the relations hold 



^'^g^' = i\/a'^[iv{x,y)] 
da{x,y) 



1 — X, 



dy 



'ix,y)]- 



(18) 
dyi, 

(19) 

(20) 
(21) 



As follows from Eq. ([19 
da{0,y)/dy 



iv{0,y) — y and therefore 
•y/7 + Q?(y). To obtain the function 



a{x,y) one has to determine the function v{x,y) from 
Eq. ini) and to insert it into Eq. (fTS)) . 

Eq. is a condition of independence of a on "con- 
stant" v{x,y). It is expressed by the relation dcr/dv = 
which is generic with the saddle condition k — ks{x,y) in 
Eq. ((T3|). We see that ks{x,y) = a[iv{x,y)] and the last 
term in Eq. (fT8|) is —F{k) if to express v through k from 
the equation k = a(w). 



VI. SOLUTION 

After a little algebra the equation reads 



a; = (t) -f iy). 



1 4- a^(iti) (w + iy) 



(22) 



Since a{y) is an even function there is a real solution 
?;(a;,0) of Eq. (P^ . Therefore there is an imaginary 
branch of (j{x, 0) as follows from Eq. 




FIG. 2: The plot of a; versus v at y — according to Eq. (122] 
at ao < a. See Sec. IVII 



We specify now the particular form 



"(y) = aoexp 



(23) 



One can easily plot x as a function of w at y = 0. At 
QfQ = 0.03, 7 = 0.2, and a = 2 the plot is shown in 
Fig. [21 It determines the real function v{x, 0). The curve 
in Fig. [2] has two special points marked as ci and C2. 
The variable x cannot be extended above the point ci 
and below the point C2 in frameworks of a real v. The 
same is valid for the action a and, therefore, for the wave 
function ^ plotted in Fig.jS^a). 

In Fig. ISIa) ci and C2 are branching points where two 
of three branches 1, 2, and 3 merge. Extensions of the 
branches are shown by the dashed curves corresponding 
to extensions of v(x, 0) to the complex plane above ci and 
below C2 in Fig. |2l Near the point C2 one can estimate 



icr{x, 0) - io-(xc2 1 0) - {xc^ 



x)\/a'^{ivc2) + 1- Xc2 
(24) 



Analogous relation holds for the point ci . 

At the points a and b the function tp{x,0) reaches ex- 
trema. According to Eq. (|20p . at the point b 

Xb^l + a^iivb), (25) 

where Vb is determined by the largest root of Eq. ([221) 



Vb = 2y^l + a'^{ivb)\/"f + a'^{ivb). 



(26) 



Equations analogous to and ^2E\i also hold for the 
point a when one should take the smaller real root of 
Eq. (|26l) . The points a and b are marked in Fig. [2l 

What happens to ip{x, y) at a finite yl To answer this 
question we consider a classical trajectory in the potential 
{—x) relating to the total energy 7 — 1 



X = Xb + 



y 



A{xb - 1 + 7) ' 



(27) 



As follows from the equation (|^ . v{x,y) ~ Vb \i x 
and y are connected by the relation ([27]). According to 
Eqs. (pOI and (HH), imaginary part of cr{x,y) is a con- 
stant on the classical trajectory (1771) . This means that 
the extreme point of the branch 3 in Fig. [3] moves in the 
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FIG. 3: Branches of the wave function at ao < a. (a) At y — 
there are two branching points ci and C2. (b) The case of a 
large y. The branches 1 and 2 at a small x are generic with 
conventional WKB ones in a one-dimensional case. The top 
of the branch 3 moves as a classical particle to infinity. 



l¥(x,0)l (a) 




FIG. 4: Branches of the wave function at a < ao. (a) Hy- 
bridized branches 1-3 and 3-1 and the branch 2-2 at y = 0. 
(b) Branches 1-1, 2-2, and 3-3 touch each other at the singular 
point ci when y — A. Another singular point C2 corresponds 
to y = -A. 



plane {x,y} along the classical trajectory ([27]) where 
keeps a constant value. This statement is valid in the 
main exponential approximation used when a quantum 
mechanical smearing of packets is negligible. 

Along the part y < of the branch 3 the particle flux 
is directed towards the well and at < y it has the op- 
posite direction. So the total flux of particles, associated 
with the branch 3, in the x direction is zero but in the 
y direction it has the same sign as the flux along the 6 
well. 

In Fig. ^h) the branches are shown at a large y when 
the 6 well is homogeneous in the y direction. In this case 
the branches 1 and 2, at a not large x, are generic with 
conventional WKB ones in one dimension. The dashed 
part of the branches 1 and 2 corresponds to an expo- 
nentially small outgoing flux away from the barrier. The 
branch 3 in Fig. ^h) is located at a large x corresponding 
to a far point of the classical trajectory (ITT)) . 

The behavior of branches of the wave function in Fig. [3] 
occurs when the width a of the 6 well profile is not too 
small. Under reduction of a a branch behavior becomes 
qualitatively different. Namely, the region between the 
points ci and C2 shrinks and becomes zero at a; = xq 
and V = vq when a = qq. Note that the points a and b 
remain well separated. For ao and 7 used in this section, 
ao — 1.72, Xq ~ 1.07, and vq ~ 1.25. When two extrema 
in Fig. [2] are about to coincide Eq. ((22l) reads 

x-xo+0.2iiy = 0.96{v-vof-0.90{a-ao){v-vo). (28) 

The singular part of a solution of the cubic equation ([^5]) 
can be evaluated as 

{v - vo) - ^/x-xo + 0.24i{y-A), (29) 

where 

A~ 1.43(ao-a)3/2. (30) 



According to Eq. (|^D|) . the part contributes to 

da/dx. Therefore at ao < a 

if7(a;,y) -icr(2;o -l-0.24|A|,0) (31) 
r^[x-xa- 0.24|A| + Q.2'iiyf''^ 

and at a < ao 

icj{x, y) - ia{xo, A) - [a: - xo + 0.24i(?/ - A)f^^ . (32) 

The singular part of the action (IHT1) corresponds to 
Eq. (|24| since Xc2 = xo + 0.24|A|. The positions of sin- 
gularities in the {x, y} plane are 

{xo ±0.24|A|, 0}, ao<a (33) 
{xo, ±0.24A}, a < uq. (34) 

The singularities ([55]) relate to Fig. [3Ka) and ones ([M]) 
relate to Fig. [3)^b). At a < ao branches 1 and 3 are 
shown in Fig. ID At y — A the branch 1 touches the 
branch 3 at the point ci resulting in the singularity (j32p . 
Analogous singular point C2 corresponds to y = —A. At 
— A < J/ < A reconnection occurs when the branches 1-3 
and 3-1 are formed. Fig. |D^a). At A < \y\ the branch 
3 detouches and its top moves as a classical particle to 
infinity analogously to Fig.[3I^b). 

VII. SEMICLASSICAL APPROXIMATION IN 
TWO DIMENSIONS 

We use above the semiclassical formalism to solve the 
Schrodinger equation p^ . This is possible when the 
semiclassical parameter B is large. In this case one can 
reduce the exact equation to the Hamilton- Jacobi 
type (fTB|) and to obtain the wave function (fTil) in the 
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Imx 



point 







Rex 



FIG. 5: Stokes lines in the complex plane of x in the case of 
a one-dimensional overbarrier reflection. Long arrows show 
an incident wave. The short arrow indicates the reflected 
wave existing at negative x only and resulted from Stokes 
phenomenon. 



exponential approximation when the action is large. To 
determine a prefactor one should go to next orders in 
the small parameter 1/B in Eq. (|15p . As a result, the 
preexponential factor is an asymptotic expansion with 
respect to 1/B. 

At the first sight, the obtained solution of the 
Hamilton-Jacobi equation provides the main (exponen- 
tial) part of the exact wave function with no problems. 
However, the prefactor may jump down to zero as a func- 
tion of coordinates. This is called Stokes phenomenon 
1281 [23. [sol considered in this section. 



A. Semiclassical approximation in one dimension 

To demonstrate the nature of Stokes phenomenon it is 
better to start with a one-dimensional static Schrodinger 
equation. We consider a known problem of overbarrier 
reflection when a particle energy E is larger than a barrier 
hight V . To be specific we choose the potential barrier in 
the form V{x) — V / cosh^ {x / a) . In classical mechanics 
overbarrier reflection is forbidden. In quantum mechan- 
ics, generally speaking, a reflected wave is not zero. The 
incident wave in WKB approximation is 



~ exp I y "^^1 



1 - 



E 



(35) 



where k — \/2mE/h is a wave vector at the infinity 

The solution (j35p holds at every x since the momentum 
\/2ra\J E — V{x) is nowhere zero. This put a reasonable 
question, how and where a reflected wave forms. To an- 
swer the question one should consider the whole complex 
plane of x. A terminal point, where E = ^^(a;), is imag- 
inary, Xc — ia arctan E/V — 1. The exact solution of 
the Schrodinger equation with the potential V{x) can be 
written in the form 

ip — exp [(^(x)] . (36) 
In the semiclassical approximation close to the terminal 



ka 

T 



2i{x-Xc)V''^ (E-V 



V 



1/4 



(37) 



The terminal point, marked in Fig. [5] as c, is an ori- 
gen of Stokes lines which are determined by the condi- 
tion Im[(^(x)] = and are shown by dashed lines. The 
point c is a branching point. Along a Stokes line an in- 
creasing (dominant) solution of the Schrodinger equation 
increases most fast and a decreasing (subdominant) solu- 
tion falls most fast compared to neighbor directions from 
the point c. According to this definition, the term "sub- 
dominant" is applicable to a branch even when it is a 
single one. 

A dominant branch is continuous when we cross a 
Stokes line. If on some Stokes line a dominant branch 
is absent, a subdominant branch is also continuous after 
crossing the Stokes line. 

The case of coexistence on a Stokes line of two types 
of branches, dominant and subdominant, is not trivial. 
Strictly speaking, it is impossible to keep an exponen- 
tially small (subdominant) branch at those x where an 
exponentially large (dominant) branch exists. As shown 
by Stokes [28| , in presence of a dominant branch a coef- 
ficient at a subdominant branch can jump after crossing 
a Stokes line. This is called Stokes phenomenon [29, 3(|. 

In our case of overbarrier reflection the incident wave 
(|35| is indicated in Fig. [5] by long arrows. The exact 
solution (|36p at real a; is a sum of the incident (dominant) 
wave and a reflected (subdominant) one 



-0(a;) = ^l)i{x) + ipRix), 



(38) 



where in the semiclassical approximation the reflected 
wave is 

^r{x) ^ d{-x) exp -Trafc (l - ^/V/E^ (39) 
exp(^-»fc^ dx.^jl-^^"^. 

The subdominant branch indicated by the short ar- 
row in Fig. [51 jumps to zero at positive x. This Stokes 
jump happens at the vertical Stokes line in Fig. [5] 

A rapid jump of a prefactor, as in Eq. occurs in a 
semiclassical solution only. An exact solution is continu- 
ous. To show that one should go beyond the semiclassical 
approximation considering a "microstructure" of a Stokes 
line. 



B. Two dimensions. 

When dimensionality is more than one, semiclassical 
approximation becomes much more complicated js^. In 
this paper we do not concern general aspects of this 
problem. We only analyze the semiclassical solution of 
Sec.[Vll 
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caustic 



FIG. 6: Classical trajectories are reflected from caustics 
marked by the dashed lines which pierce the physical plane 
{a;,Rey} at the points {xo,iA}. The trajectories and the 
caustics lie in the planes Key = ±A. The both caustics are 
tilted by the angle 6 with respect to the vertical direction. 
The point ci is the same as in Fig. |4jb). 



First, there is an analogy in properties of wave func- 
tions in one and two dimensions. In one dimension the 
wave function has the branching point (|37p in the plane 
{Rex, Imx}. In two dimensions a similar branching point 
([5^ is in the plane {x, Rey}. 

The singularity in the {x, Rey} plane can be stud- 
ied by extention of dimensionality to three dimensions, 
{x, Rey, Imy}, as shown in Fig. [6l To do that one can 
consider the action ([32|) close to the point ci in the plane 
Rey = A in Fig. [SI that is a{x,A + iif). In the plane 
Rey = A classical trajectories are reflected from the cer- 
tain lines as shown in Fig. [6l Those lines, called caustics 
[21I [soj. came to quantum mechanics from optics. Un- 
derbarrier caustics were studied in Refs. [13, [3, [l^- The 
caustics pierce the physical plane {x, Rey} at the points 
{xq, ±A}. The classical trajectories, reflected from the 
caustics, are 



x{T],h) = 



+ (1 — 7) tan^ 9 + 6tan( 
1 + tan^ e 



(40) 



{I + tan^ e){Tj ~ b) 
4(xo - 1 + 7-l-5tan6')' 

The parameter b marks different trajectories in Fig. [S] It 
is not difficult to check that the trajectory (|40)) relates to 
the energy 7 — 1. The relation dx/drj = tan0 determines 
at each b the certain point on the trajectory where it is 
tangent to the caustic. Each caustic is defined by the 
relation x — xq = j] tan 9. The caustics are straight lines 
in small vicinities of the points ci and C2 in Fig. [6l 

As in a one-dimensional case, the singular part of the 
action depends on a distance, to the power 3/2, in the 
direction perpendicular to the caustic, that is on [[x — 
xq ) cos rj—rj sin rj\^^'^ ■ Comparing this with the expression 
([32)) . where y — A + irj, one can conclude that tan 6 — 
0.24 and therefore in Fig. [Sl^ — 13.5° for the parameters 
chosen in Sec. IVII 

In a vicinity of each point of the caustic the semiclassi- 
cal action is singular having a branching point of the type 
The caustic is an origen of Stokes planes shown in 



caustic 




FIG. 7: The caustic (the dashed line) is an origen of Stokes 
planes which intersect the physical plane {x, Rej/} along 
Stokes lines. 



Fig. El The Stokes planes in the three dimensional space 
are analogous to Stokes lines in a plane. The Stokes 
planes intersect the physical plane {x, Rey} along the 
Stokes lines plotted in Fig. [71 

Now let us move along the branch 1-3 in Fig.m^a) from 
X = 0. This branch remains dominant when x is less 
than the intersection point with the branch 3-1. Then 
it becomes subdominant within the finite interval of x. 
According to Stokes theory [2§|, at some point within 
that interval a jump of the branch 1-3 down to zero is 
possible. Therefore there is no guarantee that at a larger 
X, when the branch 1-3 becomes dominant again, it would 
be the same as a formal continuation of 1-3 according to 
semiclassical formulas. 

To analyze that situation we consider the case when 
the parameter a is close to ag defined in Sec. IVII This 
means that the points ci and C2 in Fig. [6l are close to 
each other. At parameters chosen in Sec. IVII Eq. pSjl at 
J/ = yields 

^ = 2.88 [{v ~vo- 5f + 2{v - vo - S)S] , (41) 



where 5 = 0. 56-^0 — qq is a small parameter. The value 
V — Vo + S, at uq < a, corresponds to the singular point 
C2 in Figs. [Hand [3K a). According to Eq. ^U^ . at x close 
to Xo and v close to vq one can write at ?; = 



o ■ 

= 0.13 + 0.20(1; -i;o -(5). 

ox 



(42) 



Equation (|4T|) allows to integrate in Eq. (|42|) . The result 
at a < ao is 



ia{x, 0) — idixQ, 0) — 0.13(x — xq) 



(43) 



= 0.14 



{v-vo-i\S\f + jiv-vo-i\S\m 



The right hand side of Eq. (|43|) is a singular part of the 
action in terms of the variable v. According to Eq. (|4ip. 
a small {v — vq — 6) is proportional to a square root of the 
distance and therefore the cubic part of the expression 
has the same type of singularity as Eqs. ([HT|) and 
. Another singular part of the action can be obtained 
from Eq. (1^51) by the formal change \6\ — 1<5|. 
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FIG. 8: Stokes lines of the action a{x,0) which is considered 
as a function of the complex variable v at a < ao- 



Belovif we consider the case of a < oq. The singular 
part of the action, the right hand side of Eq. (HS)) . is a 
function of the complex variable v. This part becomes 
zero at the point v — vo + i\6\ w^hich is an origen of Stokes 
lines in the complex plane of the variable v. Those lines 
are determined by the condition 



Im 



{v-vo-m^ + j{v-vo-mf\s\ 



(44) 



and are shown in Fig. [8l 

There are three branches in the problem indicated in 
Fig. [3] or Fig. [3 The branch 1-3 in Fig. H^a) for all x 
relates to real v in Eq. ([22| but the branches 1-3 and 2- 
2 are determined by v with finite imaginary parts. This 
means that on the line Imv = in Fig. [5] there is only one 
branch which is 1-3 of Fig. SJa). We see that different 
branches are localized at different parts of the complex 
plane of v. Hence they cannot be classified as domi- 
nant and subdominant. In contrast, in terms of conven- 
tional coordinates all branches exist at the same x where 
dominant-subdominant competition is unclear. 

Motion along the branch 1-3 in Fig. [3)^a) from a; = 
is equivalent to sliding along the line Imi; = in Fig. [H] 
from V — 0. According to Eq. (|43p . in a small vicinity 
of vq the single (at real v) branch becomes subdominant. 
Obviously, it is continuous after crossing the Stokes line 
at the point vq in Fig. [H 

One can conclude from here that despite the branch 
1-3 in Fig. Ufa) becomes subdominant within the cer- 
tain domain of x there are no Stokes jumps there and 
the branch 1-3 goes continuously through that domain 
according to the semiclassical formulas obtained. 

In vicinities of the points a and b in Figs. [3] and H] the 
imaginary part of the action is quadratic with respect to 
(x — Xa,b)- In a one-dimensional case an action also can 
be quadratic when V{x) — E is quadratic close to some 
point X. This point is also an origen of Stokes lines in 
the one-dimensional case (29j . 

In two dimensions the extremum of Imcr at the point 
b is continued to the whole classical trajectory (|?7|) . The 
component of the vector Vcr, parallel to the trajectory, is 
continuous if to cross it and is determined by Eqs. (I^Dl) 



and ((2T|) . Therefore the perpendicular component of Vcr, 
as a part of the whole vector Vcr, is also continuous in 
crossing the trajectory. The same is valid for the point 
a. It follows that in vicinities of the points a and b the 
branches are determined by the semiclassical solution and 
look as ones shown in Figs. [3] and 21 

So we proved that the semiclassical solution of Sees. |V] 
and IVII determines the exact wave function with expo- 
nential accuracy. 



VIII. CLASSICAL TRAJECTORIES IN 
IMAGINARY TIME 

A particle penetrates through the barrier from the S 
well to the outer region where it moves along the clas- 
sical trajectory (j27|) . Accordingly, a barrier penetration 
parameter can be defined as 



w • 



ip{xb,0) 



V'(0,0) 



(45) 



In Sec.|X]we discuss a meaning of the parameter w. In the 
semiclassical approximation |^/;(0,y)| does not depend on 
y and \ilj{x,y)\ on the trajectory P7)) is also a constant, 
{xfc. 0)} is just a point of the trajectory (P7)) . By means 
of Eq. (|14p the parameter (|45)) , with the exponential ac- 
curacy, is 



W ' 



exp{-2SIm[cr(j;fc,0) -cr(0,0)]} 



(46) 



To find w one can substitute the solution of Sec. IVII into 
Eq. (gni). 

But there is another way. For a calculation of w it 
is not necessary to know the action in the entire {x, y} 
plane. If we know the action on a certain path x{y), 
which connects two points {xb,Q} and {0,0}, it is suffi- 
cient to determine w. A role of such path can be played 
by a classical trajectory in imaginary time since a mo- 
tion under a barrier in real time is classically forbidden. 
In the method of classical complex trajectories in imagi- 
nary time t — IT the coordinate x{t) remains real but the 
other coordinate becomes imaginary ?/(t) = iri^r). This 
type of complex coordinates was used in magnetotunnel- 
ing [13. See also Refs. [ll,|3|. 

The Hamilton- Jacobi equation (|16p . according to rules 
of classical mechanics [s^ , generates the equations of mo- 
tion 



2dT' 



= -1, 



1 c>\ 
2dT^ 



= 0, 



where time is measured in the units of HB/uq. 
boundary conditions to Eqs. (|47)) 



x{0) = Xb, ^ 
'7(0) =0, 

OT 



= 0, 



(47) 



The 



(48) 



= 2^Xb -1 + 7 
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FIG. 9: The bypass in the plane Key = is the solid curve 
along which the action is tracked. The part, connecting the 
points {0, iui,} and {a;i,,0}, relates to the classical trajectory 
in imaginary time. The caustics are shown by the dashed 
curves. 



are compatible with the energy conservation 



1 fdx 
The trajectory 

x{t) = Xb - 



1 fdrj 



r^, 77(t) = iT^Jxb - 1 



(49) 



(50) 



is a solution of Eq. (j47| . It starts at the point r = 
(the exit point from under the barrier) and terminates 
at To = yjxh when x(tq) = (the position of the (5 well). 

The terminal value 77(to) — 2\J x\,{xh — 1 + 7) coincides 
with Uf, as follows from Eqs. and It also follows 
from Eq. (|49p that the relation 



dx 

a7 



-2Vl + a2(wh), 



(51) 



holds which corresponds to the condition (|17p at the b 
well. 

The trajectory (|5D|) is shown in Fig. [5] as the curve 
connecting the points {a;b,0} and {0,Wfc}. The barrier 
penetration parameter (|46p can be written in the form 



where 



and 



w ~ exp(-Ao - Ax) , 
Ao = 2B Im [a{xb, 0) - ct(0, ivt)] 
Ai = 2Mm [cr(0, ivb) - cr(0, 0)] . 



(52) 
(53) 
(54) 



By means of the classical trajectory (|50)) one can cal- 
culate the part Ap as an integral of the Lagrangian with 
respect to time [36[ 



An = 2B 



1 rdx 

4 I a7 



dr] 

a7 



dr. 

(55) 



With the trajectory ([50]) Eq. ([55]) reads 

= ^ v/l + a2(iz;b) [1 - 37 - 2a2 (tVb)] . (56) 

The part Ai cannot be determined by trajectories. One 
can use directly the Hamilton-Jacobi equation at 
X — with the condition (fT7|) . One should integrate the 
term da/dy. The path, connecting the points {Q,ivb} 
and {0,0}, is the vertical solid line in Fig. (5] After the 
integration we obtain 



Ai = 2B 



(57) 



On the other hand, the same result for Ao + Ai fol- 
lows from the direct solution (fT8|) of the Hamilton-Jacobi 
equation integrating the last term by parts and using the 
relation (p6)) . 

One should note that in the case of a homogeneous 
(in the y direction) wire Eq. (j52p yields the conventional 
WKB result for tunneling through a triangular barrier 



w 



exp 



AB 



a{y) = 0. 



(58) 



The whole bypass, connecting the points {a;f,,0} and 
{0,0}, is shown in Fig. [9l It goes in the plane Rey = 
between the two caustics which are continuations of 
ones from Fig. [6l These caustics are origins of Stokes 
planes as in Fig. [T] One plane starts at the caustic ci 
and terminates at C2. Therefore the path in Fig. [9] pierces 
that Stokes plane. According to Sec. IVIH one can track 
the semiclassical solution along the path in Fig. [5] with 
no Stokes jumps. 

The conventional (no tangent momentum) underbar- 
rier trajectory in Fig. [TJa) corresponds to real coordi- 
nates. In contrast, the path in Fig. [5] lies in the complex 
plane providing a real tangent velocity dy/dt = drj/dr. 
This path can be treated as a saddle point with fluctua- 
tion around it. Those fluctuating paths in the real plane 
{x, Rey} look as a wide bundle of trajectories of the type 
shown in Fig. [Tfb). 



IX. ABOUT A NUMERICAL CALCULATION 

In principle, the Schrodinger equation ([T]) can be solved 
numerically using a discrete two-dimensional net with 
small finite steps Ax and Ay. The error in use of such 
numerical scheme is of the order of (Ax)^ ^ (Ay)^. A 
reduction of steps increases a precision of numerical cal- 
culations. 

A specificity of the problem is that one should keep 
at the same spatial domain two branches of the wave 
function, dominant and subdominant. In Fig. ^h) the 
branches 1 and 2 merge when, according to Eq. (|58p . 
I?/)! ~ exp (— 2_B/3). The subdominant branch 2 at x = 
in Fig.[3]^b) can be estimated as ~ exp (—45/3). 
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To resolve the exponentially small subdominant branch 
the precision of numerical calculations should be no lesser 
than (Ax)^ exp(-4B/3). So the number N ^ 
1/AxAy of points in the two-dimensional discrete net 
should be no smaller than exp (4i?/3). 

The parameter B has to be sufficiently large to get a 
developed semiclassical regime. On the other hand, B 
should not be too large since it results in a nonrealisti- 
cally large number of discrete points in the net. A choice 
of B for numerical calculations is a compromise between 
the two tendencies. 

According to our experience in numerical studies of 



this type of problems [27[, a semiclassical regime occurs 
when B is no smaller than approximately 25. In this 
case the number of points N should be no smaller than 
10^^ which is essentially outside of memory facilities of 
computers. This makes impossible a straightforward nu- 
merical investigation of the problem. 



X. HOW A PARTICLE PENETRATES 
THROUGH THE BARRIER 

We specify the same parameters 7 = 0.2 and = 0.03 
as in Sec. I VII One can use the direct solution of the 
Hamilton- Jacobi equation, Sec. IVH or trajectory results. 
Sec. I Vim to calculate the barrier penetration parameter 
p5)l which, at a small |afl — a|, has the form 



■w ^ exp [— 2.0i3(a_R — a)] , 
where an ~ 2.27. In dimension units Eq. ([59|) reads 



(59) 



w ~ exp 



-4.54B 1 



(60) 



where £r ~ 2.21 uq/ a. The result ((60)) holds when 
|1 — is small. As follows, when the electric field is 

close to Er the parameter w becomes not exponentially 
small. In this case in Fig.[3I^b) xi, ~ 2.0. One should note 
that tunneling rate a,i £ = Er from a homogeneous wire, 
a(y) = 0, is exponentially small according to WKB. 

The barrier penetration parameter w is rather formal. 
To physically formulate a problem of barrier penetration 
one should localize a particle in a vicinity of the 5 well 
and than study a subsequent dynamics. Suppose that the 
wave function, which is localized at the 5 well at i = 0, 
is ^p{r), where r — {x,y}. The function ^p{r) coincides 
with il){r) in Fig. [3] at x < Xa and ^p{r) = at a larger x. 

The system, being released at < = 0, starts up with 
the abrupt function ip{r) to restore the true wave func- 
tion. In the conventional case of a homogeneous wire, 
a{y) = 0, during the short time h/uf) the branches 1 and 
2, as in Fig. [3I^b), are formed. Then the process becomes 
exponentially slow when a weak leakage through the bar- 
rier is provided by an outgoing flux related to the right 
hand side of the branches 1 and 2. This process can be 
described in terms of an exponentially small imaginary 
part of the total energy. 



y 






X 







FIG. 10: The scheme of an expansion in time of the state 
which tends to restore the true wave function localized along 
the classical trajectory shown by the dashed curve. 



When the wire is not homogeneous, a{y) ^ 0, a sce- 
nario becomes different. The function 'ip{f), plotted in 
Fig. [3] is almost an eigenfunction with a real energy if to 
ignore an exponentially weak leakage across the barrier 
provided by the right hand side of branches 1 and 2. Any 
local artificial distortion of the eigenfunction is restored 
during a time interval inversely proportional to the bar- 
rier height. This is similar to the case of a homogeneous 
wire. The use of the abrupt function (p(r) instead of ex- 
act ip{r) is an example of such distortion. Therefore a 
conversion of ^{r) into is characterized by the time 
h/uo- Within the short (nonsemiclassical) time h/uo, due 
to the uncertainty principle, states with all energies are 
involved including overbarrier propagating waves. These 
waves provide a probability transfer from the S well to 
the outer region. 

The subsequent dynamics is schematically shown in 
Fig. [To] where the expanding in time probability density 
tends to restore V'('^) localized on the classical trajectory 
(l?7l) . This expansion occurs infinitely if the space Xb < x 
is not restricted. 

At a relatively small electric field, £ < £r, the parti- 
cle density on the trajectory is exponentially small. The 
decay rate (filling out of the trajectory in Fig. (TO]) of the 
state at the 5 well is exponentially weak and is deter- 
mined by ly, Eq. ||SD|) . We do not consider a preexponen- 
tial factor having the dimensionality of inverse time. 

At a larger electric field, £r < £, the peak of den- 
sity outside the barrier, at a; = x;,, is exponentially large 
compared to one at the S well as follows from Eq. (|60p . 
This means that the region at the S well is emptied, down 
to an exponentially small density, during the short time 
scale h/uQ. 

So one can describe the phenomenon as follows. Sup- 
pose a state in the non-homogeneous wire to be not a 
ground state. Under increase of the electric field, above 
the certain threshold £fi, the prebarrier region gets emp- 
tied fast (instant tunneling). The threshold £ii depends 
on energy of a decaying state in the quantum wire. For 
a higher energy £fi is less. 
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FIG. 11: The potential barrier with the impurity, related to 
Eq. (|6ip . is plotted for positive x. The state with a finite 
prebarrier momentum is shown by the arrow line. The state 
after the barrier is indicated as the classical trajectory re- 
flected from the barrier at the point x ~ 1, y — marked by 
the dot. 



XI. TUNNELING THROUGH A BARRIER 
WITH AN IMPURITY 



In this section we consider tunneling from a homoge- 
neous quantum wire described in Sec. HIl where one has to 
put a{y) = 0. But now the barrier is non- homogeneous 
in the direction parallel to the wire. Namely, there is an 
impurity at the barrier region sketched in Fig. 1111 In the 
dimensionless units of Sec. [W] the Schrodinger equation 
has the form 

^AlA + ^)+ [^0 (x) + u{x, y)] V; = (7 - 1)^, 



52 \dx 

where y) is a potential of the impurity and 

2 



B 



5{x)-\x\. 



(62) 



For simplicity we consider a symmetric impurity poten- 
tial u{—x, y) = u(x, y). This is equivalent to two identical 
impurities symmetrically localized around the wire. 



A. Hamilton-Jacobi approach 

We apply a semiclassical method based on the 
Hamilton-Jacobi equation to tunneling through a barrier 
with an impurity. The wave function has the form 
resulting in the Hamilton-Jacobi equation for positive x 



da 

dx 



da_ 

dy 



■ + u{x,y) =7-1. 



(63) 



The equation should be complemented by a 

boundary condition at a; = 0. The function a{x,y) is 
continuous on the line a; = and (da/dx)^ has the same 
values at x = ±0. Since the total potential is symmet- 
ric with respect to x the boundary condition at x = -1-0, 



accounting for the S function in Eq. ([6T|l . is 



da{x,y) 



dx 



(64) 



a:=0 



We consider in this section a weak impurity, u(x,y) 
1, which enables to treat it as a perturbation in the 
Hamilton-Jacobi equation. We emphasize that the possi- 
bility to consider u{x,y) as a conventional perturbation 
in the Schrodinger equation (j6ip is realized under more 
rigorous condition u{x,y) <^1/B. 

If to neglect u{x,y) the solution of the Hamilton-Jacobi 
equation is 



ao{x,y) — ky + i / dxi\/l — xi, 7 = fc 



(65) 



where the real parameter k has a meaning of a wave 
vector parallel to the wire. At a finite k the solution 
([65ll is not a ground state and corresponds to the finite 
underbarrier phase ky. This is an essential feature of the 
underbarrier state. The state with a finite momentum 
in the S well is indicated in Fig. [TT] as the arrow line. 
Since we consider an underbarrier state with a negative 
energy, it should he < 1. Strictly speaking, 7 contains 
an exponentially small imaginary part related to decay 
of the metastable state at the 5 well. We omit that small 
correction in calculations of the wave function. 

In the next order with respect to u{x, y) the solution 
can be written in the form 



(7 = (7o{x,y) + (7i{x,y), 



d(Ji{x,y) 



dx 



x=0 



(51) The correction cti satisfies the equation 



..dai 



d(7i 



2i——\/l - X + 2k—— = -u{x, y). 
ox oy 



0. (66) 



(67) 



The solution of Eq. (|67p , obeying the boundary condition 
dSSl), is 



CTi(a;,y) = 



[ (0, yi+y- 2ik^/l -x + 2ik) 

Jo 2k 



(68) 



idxi 



2y/l - Xl 



(a;i, y — 2ik^/\ — x + 2ik^/l — xi) . 



The correction cti to the action describes scattering of 
underbarrier waves by the impurity. 



B. Solution 

Now one should specify a particular form of the impu- 
rity potential u{x, y). We use the exponential form 



u{x,y) 



-wexp 



-u exp 



{x - If + y2 



(69) 



{x + lf +y 



21 
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with a small parameter a 1. The impurity potential is 
well localized. The dimensionless parameters Vq, I, and a 
can be easily expressed through corresponding physical 
ones. 

A simple analysis of Eq. (|68p shows that outside the 
barrier a maximum of \tp{x, y)\ is reached on the classical 
trajectory y = 2k\/x — 1. In the semiclassical approxi- 
mation used \^{x,y)\ is a constant along the trajectory. 
It smears out if to account for quantum effects described 
by the last terms in Eq. (|15p . That solution, localized at 
the classical trajectory, is symbolically shown in Fig. [11] 
by the arrow curve. 

The integration in Eq. ([68]) is not difficult. Under the 
condition 



I <2k^ <2 



(70) 



the modulus of the wave function, at 1 < x and close to 
the classical trajectory y ~ 2k\Jx — 1, has the form 



\'^{x,y)\ 



exp 



-IBj-i 



(x- 1)1/4 "^""P ] 8fc2(2F~0 



exp 



(71) 



iy - 2kVx - 1 
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COS 



4k 

'2 



— (y - 2fcV^^) 



At u = (no impurity) Eq. (|7ip turns into a conventional 
WKB expression. Eq. ([7T]) describes the state outside 
the barrier indicated by the arrow curve in Fig. [4] This 
state is generic with one at the right hand side part of 
Fig. [Tfb). The state ([7T]) . driven by underbarrier mech- 
anisms, essentially differs from a state outside the bar- 
rier without the impurity (incident and reflected waves) 
which is hardly influenced by the d well. 

Far from the impurity, \y\ — > oo, and away from 
the classical trajectory the wave function coincides with 
the conventional WKB form which is the first factor in 
Eq. ([7T]) . In this case there is an outgoing wave only out- 
side the barrier. It provides an exponentially small imag- 
inary part of the energy which is neglected in Eq. f7T]) . If 
to ignore the exponentially weak leakage across the bar- 
rier the wave function (|7ip can be treated as an eigen- 
function. 

Above we did not account for a small correction 5^ to 
the eigenvalue 7 — 1 (|65p due to the impurity potential. 
This real correction would result in the additional part 
yS^ /2k in (Ji{x, y) which does not influence the modulus 
of the wave function ([71]l . 

The applicability conditions of the result ([7T]) are 

uexp( — — — j < 1, exp ( — — — j<B. (72) 

We remind that i? is a large parameter. The first inequal- 
ity (|72p follows from the perturbation condition with re- 
spect to the impurity potential. The second condition 
(17^ is semiclassical one when the last part in Eq. is 
less than u{x,y). 
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FIG. 12: The processes of scattering by the impurity placed 
close to the position a; = 0. The eigenstate in the total poten- 
tial Vo{x) + u{x, y) (the barrier plus the impurity) corresponds 
to the underbarrier eigenvalue of the energy 7 — 1. The prop- 
agating waves (eigenstates of Vo{x)) are a part of the total 
eigenstate. The dots mark impurity scattering. 



C. Results 

One can draw two conclusions on the basis of equation 
([7T]) : (i) the effective amplitude of the impurity potential 
Mexp[(4fc^ — P)/a^] is exponentially enhanced compared 
to u since fc ~ Z ^ 1 and a is small and (ii) a scenario of 
barrier penetration is of the type as in Fig.[T]^b). 

The underbarrier exponential enhancement is generic 
with one occurring in tunneling through nonstationary 
one-dimensional barriers (25j where there is an interfer- 
ence of various paths generated in different moments of 
time. 

In Fig. [TT] the impurity position / is before the exit 
point. According to the conditions (|70p . the impurity 
can be placed even after the exit point, 1 < L 



XII. THE NATURE OF THE PHENOMENON 

In this section general arguments are formulated which 
allow to analyze the phenomenon by "hand waving" . To 
be specific we consider tunneling from a homogeneous 
quantum wire through the barrier with an impurity (|6ip . 

In absence of the impurity {u(x,y) = 0) motions in x 
and y directions are independent. In this case tunneling 
in the x direction occurs according to the WKB one- 
dimensional scenario. 

The eigenstate with the underbarrier energy 7 — 1 in 
the total potential Vo{x) + u{x, y) can be considered as a 
superposition of eigenstates of the potential Vq (x) : (i) the 
underbarrier WKB type state and (ii) propagating waves 
with overbarrier energies. In the first order of a pertur- 
bation theory with respect to u{x, y) the processes (a) in 
Fig- [H] participate in formation of the exact underbarrier 
state. Subsequent scattering processes, of next orders of 
the perturbation theory, modify the propagating waves 
as shown in Fig. [T^ by the paths (b) . The term "eigen- 
state" is used since one can ignore an exponentially small 
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leakage across the barrier. 

A principal question is that how strong is a contribu- 
tion of the propagating waves to the total wave function. 
First of all, if u{x, y) is not small, each partial propagat- 
ing wave is also not small and does not decay exponen- 
tially with distance. But there is an opposite tendency. 
When the waves are distributed smoothly in energy and 
directions the mutual interference reduces their contribu- 
tion to the wave function down to an exponentially small 
value. 

However, a completely different scenario can be real- 
ized when a tangent momentum in the prebarrier region 
is not zero as in Fig. [Ub). In our case this corresponds 
to a finite momentum py in the wire (x = 0). The finite 
tangent momentum may result in a new phenomenon in 
the distribution of propagating waves. A classical tra- 
jectory in the potential (— x) (at a not small positive x 
where the impurity potential is weak) can be generally 
written as 



+ 1 - 7 + 



(73) 



Here ib is a reflection point in y, py is a tangent momen- 
tum, and the underbarrier energy 7 — 1 is negative. The 
trajectory (|^ is a particular case of the form ([75]) . 

A type of the trajectory ([75)1 strongly depends on py. 
At zero tangent momentum the trajectory (j73p is degen- 
erated into the line y — when 6 = 0. This is the main 
underbarrier path of the same type as in Fig. [Ija) . At a 
flnite Py the curves (j73p become two-dimensional. They 
are analogous to rays in geometrical optics. When the 
optical rays are not parallel they are reflected by cer- 
tain curves which are caustics [2l|. The same caustic 
phenomenon in the space {x, Imy} occurs with the tra- 
jectories ([75)) as follows from the results of Sees. I VII and 

Eni 

From the standpoint of the three-dimensional space 
{a;, Hey, Imy}, the caustic curve pierces the physical two- 
dimensional plane at the point {xq, A} as shown in Fig.[6l 
In a vicinity of this point a distribution of the propagat- 
ing waves becomes not smooth (as close to a caustic in 
optics) and their interference cannot now lead to a mu- 
tual compensation in contrast to a smooth distribution. 
As a result, a significant fraction of the prebarrier den- 
sity is carried away by propagating waves. Therefore, a 
barrier penetration is determined in that case not by a 
conventional underbarrier mechanism, as in Fig.[TI^a), but 
by interference of the propagating waves, as in Fig.[TI^b)). 

We outline by general arguments the nature of unde- 
barrier interference. To draw more exact conclusions we 
use the rigorous method of summation of various prop- 
agating waves in Fig. [T^ This method is generic with 
a saddle formalism which allows to collect rapidly os- 
cillating functions. In our case this is the semiclassical 
approach based on the equation of Hamilton- Jacobi used 
in Sec. IVl 



XIII. TWO ENTANGLED PARTICLES 

A process of tunneling in two dimensions is sketched 
in Fig. [TJ This tunneling scenario can also be realized in 
a two-dimensional potential barrier of the type 



V{x,y) 



V 



V 



cosh (x/a) cosh (y/a)' 



(74) 



when the particle energy E is less than V. Suppose tun- 
neling to occur symmetrically, that is in the direction 
X = y in the plane {x,y}. A tangent momentum is in 
the direction x = —y, which is perpendicular to tunnel- 
ing. One can show (we will present calculation elsewhere) 
that in this case a phenomenon of underbarrier interfer- 
ence also takes place. The phenomenon is generic with 
one investigated above and resulting in the enhanced tun- 
neling. 

The problem of a particle in the two-dimensional po- 
tential (j74p can be considered from another standpoint. 
Namely, there are two different particles described by 
the coordinates x and y and moving in identical one- 
dimensional potentials. In classical mechanics the two 
particle are completely independent and can be spatially 
separated by an arbitrary long distance. 

In quantum mechanics everything depends on a type 
of the wave function. The wave function of two particles 
ijj{x,y) — '4'in(x)ip2n('y) dcscribcs independent motions 
through each separate barrier. Suppose now that the 
two particles are in the entangled state 



i1pln{x)')jj2n{y) 



(75) 



related to a momentum of the total system in the direc- 
tion X = —y and to some classical turning point along 
the line x = y. The wave function (j75p is analogous to 
one for a single particle in the two-dimensional potential 
(|74)) with the total energy E. Therefore, the two parti- 
cles can penetrate the barrier instantaneously and with 
a not small probability despite they are separated by an 
arbitrary distance. So by means of the second particle 
the conventional WKB mechanism can go over into one 
with the underbarrier interference. 

This scheme provides some sort of quantum commu- 
nication when one can induce a penetration through a 
barrier of the "controlled" particle which is distant from 
the "controlling" one. A role of particles with barriers 
can be played by two Josephson junctions 39]. With no 
interference effects tunneling in Josephson junctions was 
initially observed in Ref. [40]. We do not discuss here a 
way to create the properly entangled state of two junc- 
tions. 



XIV. DISCUSSIONS 

The counterintuitive phenomenon of easy penetration 
through barriers which are nontransparent, according to 
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WKB theory, is proposed. A magnetic field is zero, the 
potential barrier is static and non-homogeneous in the 
direction perpendicular to tunneling. 

When the applied electric field is less than the certain 
value decay of the prebarrier state is exponentially 
weak and, in principle, similar to a conventional WKB 
decay. 

At a larger field, £r < £, a scenario of barrier pene- 
tration becomes different. An initially created prebarrier 
state will be transferred fast (during the time of an in- 
verse barrier height) outside the barrier. 

An essential point of the phenomenon is a generation 
of propagating waves under the barrier when, prior to 
tunneling, a particle has a momentum perpendicular to 
a tunneling direction (a tangent momentum) . Tunneling 
from the ground state is of a conventional WKB nature 
since a tangent momentum is zero in this case. In pres- 
ence of a finite tangent momentum the propagating waves 
do not cancel each other completely by interference and 
the surviving part provides a not small output from un- 
der the barrier. 

The state outside the barrier, described in Sec. lVIi is a 
static packet. The top of the packet follows the classical 
trajectory ([27| where it keeps a constant amplitude. If to 
account for quantum effects, beyond the Hamilton- Jacobi 
approach, its amplitude reduces at a larger \y\. The mo- 
mentum, associated with the packet, is gained from the 
quantum wire. The momentum is directed towards the 
barrier at y < and it is opposite at < y. The packet 
is generated at the points {a;o,±A} where the caustics 
pierce the physical plane {x,y}, Fig. IHb). At A < \y\ 
the packet outside the barrier and the initial branch 1-1 
are disconnected. At \y\ < A the initial branch softly 
undergoes into the packet, the curve 1-3 in Fig.[4l^b). 

An exact analytical solution is not obtained. However, 
on the basis of the Stokes theory it is proved that the 
semiclassical exponent is a main part of the exact wave 
function in the problem. A prefactor does not unexpect- 
edly jump and, therefore, it plays a secondary role. 

Suppose the underbarrier state to be artificially cut 
off at the initial moment so that there is no the outer 
packet. After a release the dynamical state is developed, 
when transitions in the entire spectrum occur. The typ- 
ical energy scale involved is the barrier height uq and 
therefore the time scale of restoring of the branch 1-3 
is of the order of H/uq. This is a short (nonsemiclassi- 
cal) time. Under increase of the electric field, above the 
certain threshold £r, the prebarrier region gets emptied 
fast (instant tunneling). Se. depends on energy of a de- 
caying state in the quantum wire. For a higher energy 
the threshold Er is less. 

As we see in Sec. IXI[ the scattering of underbar- 
rier waves by the impurity also results in the wave 
packet propagating outside the barrier. The both situa- 
tions, described in Sees. IVII and IXIl are generic and can 
be unified as scattering of underbarrier waves by non- 
homogeneities. The phenomenon of underbarrier inter- 
ference provides a different aspect in study of tunneling 



through barriers with impurities [3l|, Is^, [H, [s^l since an 
individual impurity may strongly increase the tunneling 
rate. 

Tunneling through a barrier with an impurity can be 
unusual even for zero tangent momentum of a particle 
at a prebarrier region. But in this case the impurity 
should be dynamic. A typical example is alpha decay of 
a nucleus assisted by an incident proton ^] . The moving 
proton plays a role of a dynamic impurity. Caustics are 
formed in the joint space of alpha particle and proton 
coordinates. Due to interference in the {ck,p} system, 
moving protons can increase the alpha decay rate making 
it not exponentially small. 

Tunneling across a one-dimensional barrier with a 
nonstationary slope also can be strongly enhanced [25j . 
Quanta absorption and emission result in an underbar- 
rier phase analogous to one produced by the prebarrier 
momentum py in two dimensions. This leads to a phe- 
nomenon, analogous to caustics in Fig. [51 which prevents 
a mutual cancellation of the propagating states. 

The underbarrier interference is also a feature of tun- 
neling in multi-dimensional systems. An elastic string in 
a washboard potential can tunnel from some valley (the 
initial valley) to another one (the basic valley) through a 
potential barrier [ul, [H, [H, llllil E^l . If the string has 
a momentum along the initial valley and the string or the 
washboard is not homogeneous in that direction, an en- 
hancement of tunneling, caused by the interference, may 
occur. Calculations will be published elsewhere. As a re- 
sult of tunneling, the certain part of the string appears in 
the basic valley and extends along its direction. In pres- 
ence of a friction the part, coming to the basic valley, 
extends with a constant velocity. From the standpoint of 
the basic valley, this phenomenon looks as a violation of 
energy conservation if the basic valley is thought to be a 
single one in the general potential. 

The enhanced penetration through a classical barrier, 
caused by the underbarrier interference, can be observed 
in various natural and artificial tunneling systems. Below 
some of them are mentioned. (1) Tunneling from a wire 
or a film. Their non-homogeneities should satisfy the cer- 
tain not very rigorous conditions of the type described 
in Sec. IVII (2) Tunneling from a homogeneous wire or 
film through a barrier with impurities. (3) A system of 
two coupled Josephson junctions can manifest a not weak 
tunneling across a non-transparent barrier as a result of 
the underbarrier interference. (4) Tunneling in two dis- 
tant systems which are properly entangled (Sec. IXIII|1 . 

(5) Ionization of a molecule in an electric field can ex- 
hibit the phenomenon of interference if the molecule ro- 
tates around an axis perpendicular to the electric field. 

(6) The proton assisted alpha decay of nuclei is also an 
example of an experimental investigation. 

We studied above a pure Hamiltonian system. An in- 
terference in this case differs from one in presence of fric- 
tion, for example, phonons [12, IS]- A role of friction is 
worth to be studied. 
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XV. CONCLUSIONS 

Quantum tunneling through a two-dimensional static 
barrier becomes unusual when a momentum of an elec- 
tron has a tangent component with respect to a border 
of the prebarrier region. If the barrier is not flat a frac- 
tion of the electron state is waves propagating away from 
the barrier. When the tangent momentum is zero a mu- 
tual interference of the waves results in an exponentially 
small outgoing flux. The finite tangent momentum de- 
stroys the interference due to formation of caustics by the 
waves. As a result, a significant fraction of the prebarrier 



density is carried away from the barrier providing a not 
exponentially small penetration even through an almost 
classical barrier. The total electron energy is well below 
the barrier. 
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